Critical behavior of roughening transitions in parity-conserving growth processes 
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We investigate a class of parity-conserving solid-on-solid models which describe the growth of an 
interface by the deposition and evaporation of dimers. As a key feature of the models, evaporation of 
dimers takes place only at the edges of terraces, leading to a roughening transition between a smooth 
and a rough phase. We consider several variants of growth models in order to identify universal and 
nonuniversal properties. Moreover, a parity-conserving polynuclear growth model is proposed. All 
variants display the same type of universal critical behavior at the roughening transition. Because 
of parity-conservation, the critical behavior at the first few layers can be explained in terms of 
unidirectionally coupled branching annihilating random walks with even number of offspring. 



PACS numbers; 64.60.Ak, 05.40. -Hj, 82.20.-w. 
I. INTRODUCTION 

In the present work we continue the investigation of 
a certain class of nonequilibrium models for interfacial 
growth by adsorption and desorption Q. The models 
may be used to describe a layer-by-layer growth process 
of a d-dimensional surface by deposition and evapora- 
tion of dimers which are aligned with the surface. Upon 
adsorption the dimers dissociate into two atoms at neigh- 
boring lattice sites. It is assumed that the atoms cannot 
diffuse on the surface. Furthermore, atoms cannot evap- 
orate from the interior of completed layers. Only pairs 
of neighboring atoms at the edges of terraces are able to 
form a dimer and evaporate back into the gas phase. 

The most interesting property of this class of growth 
processes is the emergence of a roughening transition 
from a smooth to a rough phase at a certain critical ra- 
tio of the adsorption and desorption rates |^]. In con- 
trast to equilibrium growth models and nonequilibrium 
growth processes described by the Kardar-Parisi- Zhang 
(KPZ) equation where roughening transitions take 
place only in d > 2 dimensions, the models discussed 
in the present work exhibit a robust roughening transi- 
tion even in one spatial dimension. As will be shown 
below, the scaling properties of the interface at the tran- 
sition differ significantly from the usual scaling laws for 
roughening interfaces, i.e., the models exhibit anomalous 
roughening properties. 

The motivation to consider growth processes of dimers 
originates in recent studies of similar interface models 
for adsorption and desorption of monomers In these 
models individual atoms are adsorbed with probability p 
on each lattice site, whereas atoms evaporate with prob- 
ability 1—p solely at the edges of terraces. It was shown 
that monomer models of this kind exhibit a roughening 
transition which is closely related to the universality class 
of directed percolation (DP) A more detailed analysis 
revealed that the critical behavior of the first few layers 
can be explained in terms of unidirectionally coupled DP 
processes H . It turned out that such hierarchies of cou- 



pled DP processes not only describe the monomer models 
of Ref. but also various other DP-related growth pro- 
cesses, including polynuclear growth (PNG) models [Q, 
certain models for fungal growth Q , and a recently intro- 
duced model designed for real-space renormalization 
Thus, the concept of 'coupled DP' characterizes a whole 
universality class of roughening transitions ]lO[ | . 

DP itself is the generic universality class for phase tran- 
sitions into trapped (absorbing) states and is known to 
be extremely robust with respect to the choice of the 
dynamic rules. The DP conjecture |l^] states that in a 
random process with short-range interactions any tran- 
sition from a fluctuating active phase into a single ab- 
sorbing state should belong to the DP class, provided 
that the dynamics is characterized by a single-component 
order parameter without additional symmetries. Non- 
DP critical behavior is expected in systems where one 
of these requirements is violated. An important ex- 
ample is the so-called parity-conserving (PC) universal- 
ity class for phase transitions into absorbing states in 
which the parity of the particle number is conserved. In 
one dimension, this conservation law can also be inter- 
preted as a Z2-symmetry between two different absorbing 
states. The PC class is represented most prominently 
by branching annihilating random walks with two off- 
spring (BAW2) p^-p^. Other examples include nonequi- 
librium kinetic Ising models [p^, interacting monomer- 
dimer models p6| , a s well as models with two symmetric 
absorbing states |0]. It is therefore near at hand to in- 
vestigate the question how the physical properties change 
if the DP mechanism of the growth models in Ref. Q is 
replaced by a parity-conserving dynamics. To this end 
we modify the dynamic rules of these models by using 
dimers instead of monomers which adsorb with proba- 
bility p and desorb at the edges of terraces (including 
solitary dimers) with probability 1 — p. As dimers con- 
sist of two atoms, the number of particles at each height 
level is conserved modulo 2. As shown in Ref. ||l[], the 
model exhibits a robust roughening transition which is 
related to the PC class. 
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FIG. 1. Variant A in d = 1 dimension: Dimers are ad- 
sorbed with probability p and desorbed at the edges of ter- 
races with probability 1 — p. Evaporation from the middle of 
plateaus is not allowed. 



In the following we present a detailed numerical anal- 
ysis of the critical behavior of parity-conserving growth 
processes. One of our aims is to clearly identify univer- 
sal and nonuniversal properties. To this end we consider 
four variants of the model, a restricted version with ran- 
dom sequential updates, an unrestricted version as well 
as the corresponding variants with parallel updates. It 
turns out that all variants exhibit the same type of criti- 
cal behavior at the roughening transition. In particular, 
the critical behavior at low-lying levels can successfully 
be described in terms of unidirectionally coupled PC pro- 
cesses, generalizing the concept of 'coupled DP'. Above 
the transition, however, restricted and unrestricted vari- 
ants display different properties. 

The paper is organized as follows. In Sect. ^ we de- 
fine four variants of the dimer growth model. Their 
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phenomenological properties are described in Sect. HI. 
Sects. 1^ and ^ investigate the critical properties of the 
interface width and the densities of exposed sites at the 
first few layers. We also discuss aspects of spontaneous 
symmetry breaking and unusual scaling properties for 
random initial conditions. The critical behavior at the 
roughening transition can partly be explained in terms 
of unidirectionally coupled PC processes for which we 
propose a field-the oreti c formulation in Sect. VI. As will 
be shown in Sect. VII, it is even possible to construct 



a polynuclear growth model which belongs to the same 
unive rsalit y class. Our conclusions are summarized in 

Sect. Ivm. 



II. MODELS DEFINITION 

The class of models may be introduced in terms of a 
d-dimensional interface which evolves by adsorption and 
desorption of dimers. As a key feature desorption may 
only take place at the edges of a plateau, i.e. at sites 
which have at least one neighbor at a lower height. The 
dimer growth model is defined on a d-dimensional square 
lattice with L'^ sites and periodic boundary conditions. 
Each site i is associated with an integer height variable 



FIG. 2. Variants C and D: Cyclic parallel updates on four 
different sublattices in a one spatial dimension. 



hi S Z. We consider four variants A,B,C,D of the model 
which differ by their dynamic rules. 

Variant A is a restricted solid-on-solid (RSOS) model 
evolving by random sequential updates. For each at- 
tempted update a pair of adjacent sites i and j is selected 
at random. If the heights hi and hj are equal a dimer is 
adsorbed with probability p 
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or desorbed with probability 1 — p 
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min hk < hi 

ke<i,j> 



(1) 



(2) 



where k runs over the nearest neighbors of sites i and j. 
An attempted update is rejected if it violated the RSOS 
constraint 



\h,^hj\ < 1 



(3) 



i.e., the heights at neighboring sites may differ by at most 
one step. A Monte-Carlo sweep consisting of L'^ local 
update attempts corresponds to a time increment t — > 
t + 1. For d = 1 the dynamic rules are shown in Fig. ^. 
Notice that the rules are translationally invariant in time, 
space as well as in height direction. Thus, the layer at 
h — has no particular physical meaning, although we 
will often use a flat interface at zero height as initial 
condition. 

Variant B is an unrestricted solid-on-solid (SOS) model 
which is defined by the same rules without restriction (|^). 
Although unrestricted growth is less realistic, it exhibits 
essentially the same type of critical behavior at the 
roughening transition, supporting the claim of univer- 
sality. 

Variants C and D are counterparts of A and B em- 
ploying parallel updates. More precisely, their lattice is 
divided into several sublattices in a way that synchronous 
updates on each sublattice according to the rules (|l|) 
and (H) do not overlap. In d = 1 dimension at least 
three sublattices are needed. However, for technical rea- 
sons it is more convenient to work with four different 
sublattices, as illustrated in Fig. |^. We implement this 
update scheme on a parallel computer with 24000 string 
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FIG. 3. Typical interface configurations of the restricted 
dimer model (variant A) for various values of p (see text). 



processors. Associating an individual processor with each 
lattice site we are able to perform efficient simulations. 
Further details on massive parallel computing are given 
in Ref. ITl. 



III. PHENOMENOLOGY OF THE 
INTERFACE DYNAMICS 

Although the dimer models are defined in arbitrary 
spatial dimensions, this work is restricted to the one- 
dimensional case d — 1. Clearly, the morphology of the 
interface depends on the growth ratep. Let us first con- 
sider the restricted variants (see Fig. ||). If p is very small, 
only a few dimers are adsorbed at the surface, staying 
there for a short time before they evaporate back into the 
gas phase. Thus, the interface is anchored to the actual 
bottom layer and does not propagate. As p increases, a 
growing number of dimers covers the surface and large 
islands of several layers stacked on top of each other are 
formed. Approaching a certain critical threshold Pc the 
mean size of the islands diverges and the interface evolves 
into a rough state. 

Above Pc one may expect the interface to detach from 
the bottom layer in the same way as the interface of 
monomer models starts to propagate in the supercritical 
phase. However, since dimers are adsorbed at neighbor- 
ing lattice sites, solitary unoccupied sites may emerge. 
These pinning centers prevent the interface from mov- 
ing and lead to the formation of 'droplets' (see Fig. ||). 
Due to interface fluctuations, the pinning centers can 
slowly diffuse to the left and to the right. When two 
of them meet at the same place, they annihilate and a 
larger droplet is formed. Thus, although the interface 
remains pinned, its roughness increases continuously. 

As realized in a second transition takes place at 
p — 0.5 where the width grows most rapidly. In this 
case the droplets reach an almost triangular shape with 
unit slope, i.e., the surface becomes faceted. The tilted 



FIG. 4. Typical interface configurations of the unre- 
stricted dimer model (variant B) for various values of p. At 
the transition p = 0.5 large spikes of stacked dimers are 
formed. For p > 0.5 these spikes are biased to grow with 
constant velocity. 



surface of the droplets fluctuates predominantly by ran- 
domly moving triplets of sites at equal height along 'stair- 
cases' of unit slope. Inspecting the dynamic rules it is 
easy to verify that these 'landings' move upwards with 
probability p and downwards with probability 1—p. This 
explains why the faceting transition takes place exactly 
at p = 0.5. For p > 0.5 the fluctuations are confined to 
an exponentially small region at the top of the droplets. 
Therefore, the faceted interface coarsens on a logarithmic 
time scale. A pathological situation emerges for p = 1 
where evaporation of dimers is forbidden. As the pinning 
centers, once formed, cannot diffuse the interface quickly 
evolves into a frozen configuration. 

To summarize, the restricted variants A and C display 
three different phases, a smooth phase p < Pc, a rough 
phase Pc < p < 0.5, and a faceted phase p > 0.5. The 
phase structure of the unrestricted variants B and D is 
very similar (see Fig. ^). They too exhibit a roughening 
transition at a certain critical threshold Pc- However, the 
rough phase and the transition at p = 0.5 are different 
in character due to the formation of spikes, as will be 
described below. 



IV. CRITICAL BEHAVIOR OF THE 
INTERFACE WIDTH 



The morphology of a growing interface is usually charac- 
terized by its width 



h,{t) 



1/2 



(4) 



In order to investigate the scaling properties of the width, 
we perform Monte-Carlo simulations starting from a flat 
interface ft.i(0) = and get the following results: 

a. The smooth phase p < pc- Fig. || shows the 
temporal evolution of [t) for different values oip <pc 
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FIG. 5. Left: The squared width W'^(t) measured in 
variant A for various values of p < pc- Right: Data 
collapse according to Eq. visualizing the scaling func- 
tion F . Throughout the whole paper t is measured in units 
of Monte-Carlo sweeps. 



measured in model A with L — 10000 sites. As can be 
seen, the width first increases as Wif) ~ Vint until it 
saturates at some constant value. The initial logarithmic 
increase suggests the scaling form 



W'^{t,t) ~ a In tFite^i^) 



(5) 



where a is an amplitude factor and e = |p — Pc| denotes 
the distance from criticality. The exponent z^y describes 
the singular behavior of the temporal correlation length 
^11 ^ £""11 close to the roughening transition. F is a 
universal scaling function with the asymptotic behavior 
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1 const if 





oo 



(6) 



In order to estimate a and we adjust these quantities 
in a way that the curves for exp{W^/a)/t versus te''^ 
collapse onto a single one. Using the estimates Pc — 
0.317(1), = 2.4(4), and a = 0.17(1) we obtain a fairly 
convincing data collapse (see right hand graph of Fig. ||) . 
However, we observe considerable deviations for larger 
values of e. These deviations may indicate corrections of 
the scaling form ( p|) in the off-critical regime and will be 
analyzed in Sect. VE. Similar results are obtained for 
variants B,C, and D. 

b. The roughening transition at p — pc'. At 
the roughening transition the interface width of an infi- 
nite system grows logarithmically as W ^ Vint. A sim- 
ilar logarithmic behavior has also been observed at the 
roughening transition of monomer models Q]. The log- 
arithmic time dependence suggests the finite-size scaling 
form ioll 



W^{L,t) ~ a In tG{t/L^) 



(7) 



where z denotes the dynamic exponent. G is a univer- 
sal scaling function with the same asymptotic behavior 
as in Eq. (^). Thus, the interface of finite systems is 
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FIG. 6. Finite-size scaling of the interface width W{L,t) 
for variants A-D. The graphs show data collapses according 
to the scaling form (^), visualizing the scaling function G. 

expected to saturate at a constant value W{L) ^ ^/\nL. 
Our numerical results are summarized in Fig. ^ Plotting 
exp(l/F^/a)/t against t/L^ in the time interval 10 < t < 
10^ for system sizes L = 32, 64, . . . , 4096, we determine 
a and z by data collapse. For all variants we obtain ac- 
curate data collapses. The estimates for z coincide in all 
cases, indicating universal properties of the roughening 
transition independent of the RSOS constraint and the 
type of updates (see Table |). In fact, as we will argue in 
the following Section, z ~ 1.75 is the dynamic exponent 
of the PC universality class. The amplitude a, however, 
is nonuniversal. A very accurate data collapse is obtained 
by simulations of variant C on a parallel computer and 
highly supports the validity of the scaling form (|^). It 
would be interesting to confirm this scaling form by di- 
rect diagonalization of the transfer matrix plf . 

c. The rough phase pc < p < 0.5: In the rough 
phase the formation of pinning centers prevents the in- 
terface from propagating. By diffusion and pairwise an- 
nihilation of the pinning centers the width increases very 
slowly. Since the width displays a rather inconclusive 
scaling behavior, it was impossible conjecture an appro- 
priate scaling form ||2^ . It seems that the width initially 
increases algebraically until it slowly crosses over to a 
logarithmic increase W{t) ~ Valni. The amplitude a 
and the crossover time grow with p and diverge at the 
transition at p = 0.5. The crossover time provides a 
typical time scale of the dynamics in the rough phase. 
Apart from that the roughening transition is associated 
with another temporal correlation length ^|| = (p — Pc)"^" ■ 
Thus, the rough phase is characterized by a complicated 
interplay of at least two different time scales. 

d. The transition at p = 0.5: The restricted as 
well as the unrestricted variants undergo a second phase 
transition at p = 0.5 where the width increases alge- 
braically with time as W ^ t^ . Let us first consider the 
restricted case where the transition can be interpreted 
as a faceting transition between a rough and a faceted 
phase ijl^. As the width increases algebraically, we ex- 
pect ordinary Family- Vicsek scaling |23] of the form 



W{L,t)^ L"f{t/L^/^). 



(8) 
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FIG. 7. Power-law roughening of the width in large sys- 
tems. The figure shows the squared width as a function of 
time at the transition p = 0.5 (solid lines), as well as above the 
transition for p = 0.6 (dashed line) and p = 0.7 (long-dashed 
line). The inset shows the squared local width of variant A 
for increasing box sizes Z = 2, 4, 8, 16, ... as a function of time 
(see text). 



where a and /3 are the roughening and the growth ex- 
ponents, respectively. However, finite size simulations 
(not shown here) reveal a more complex behavior. After 
an initial short-time regime we observe a long transient 
extending over two decades in time where the interface 
roughens with the exponents a ~ 0.5 and (3 ~ 0.45. Af- 
ter approximately 5000 time steps both variants A and C 
cross over to a different regime with (3 ~ 0.33. By a finite 
size data collapse we find a roughening exponent a ~ 1.2 
and a large dynamic exponent z = a 1(3 ~ 3. These 
results are consistent with the findings of Ref. |l9j. In 
Fig. the crossover appears as a slight change of the 
slope, suggesting the faceting transition e.t p — 0.5 not 
to be fully scale-invariant but characterized by a finite 
time scale. Moreover, as shown in the inset of Fig. |^, 
the faceting transition displays anomalous scaling prop- 
erties [|4j, i.e., the local width w{l,t) measured in boxes 
of sizes I = 2,4, 8, 16, . . . does not saturate after a short 
time. Even for small box sizes it continues to increase 
over several time decades until it saturates due to the 
RSOS constraint. 

Let us now consider the unrestricted variants B and D. 
As can be seen in Fig. |[ their phase transition aX p = 0.5 
is different in character. Since large spikes are formed, 
the surface roughens much faster with a growth exponent 
of /? ~ 0.5 (see Fig. 0). In fact, the interface evolves into 
configurations with large columns of dimers separated by 
pinning centers. These spikes can grow or shrink almost 
independently. Thus, the interface roughens by a purely 
diffusive mechanism as W{t) ~ As the columns are 
spatially decoupled, the width does not saturate in finite 
systems, i.e., the dynamic exponents a and z have no 



variant 


A 


B 


C 


D 


restriction 


yes 


no 


yes 


no 


updates 


random 


random 


parallel 


parallel 


Pc 


0.3167(2) 


0.292(1) 


0.3407(1) 


0.302(1) 


a 


0.172(5) 


0.23(1) 


0.162(4) 


0.19(1) 


z 


1.75(5) 


1.75(5) 


1.74(3) 


1.77(5) 


So 


0.28(2) 


0.29(2) 


0.275(10) 


0.29(2) 


Si 


0.22(2) 


0.21(2) 


0.205(15) 


0.21(2) 


S2 


0.14(2) 


0.14(3) 


0.13(2) 


0.14(2) 




0.765(10) 


0.765(10) 


0.753(3) 


no result 


a 


1.2(1) 


undefined 


1.25(5) 


undefined 


/3 


0.34(1) 


0.50(1) 


0.330(5) 


0.49(1) 



TABLE I. Numerical estimates for the four variants of the 
dimer model at the roughening transition p = pc (upper part) 
and the transition p — 0.5 (lower part). 



physical meaning. 

e. The faceted/free phase: For p > 0.5 the re- 
stricted models A and C evolve into faceted configura- 
tions (see Fig. ^). As shown in Ref. the width first 
increases algebraically until the pinning centers become 
relevant and the system crosses over to a logarithmic in- 
crease of the width. Therefore, the faceted phase may 
be considered as a rough phase. The unrestricted models 
B and D, however, evolve into spiky interface configura- 
tions. The spikes are separated and grow independently 
by deposition of dimers. Therefore, the interface width 
increases linearly with time, defining the free phase of the 
unrestricted models. 



V. CRITICAL PROPERTIES OF THE 
FIRST FEW LAYERS 

A. Relation to the PC class 

In this Section we investigate the roughening transition 
at p = Pc in more detail. In order to understand its rela- 
tion to the PC class let us consider all sites at the bottom 
layer hi = as A-particles. Adsorption and desorption 
processes correspond to certain effective reactions of the 
^-particles. For example, the adsorption of a dimer at 
the bottom layer corresponds to a pair-annihilation pro- 
cess 2A — > at rate p. Similarly, when a dimer evapo- 
rates, two A-particles are created. However, since dimers 
can only evaporate at the edges of terraces, this process 
always requires the presence of another neighboring A- 
particle, giving rise to an effective reaction A — > 3A at 
rate 1—p. These two processes compete one another and 
resemble a branching-annihilating random walk with two 
offspring (BAW2) |n2 13 1 which is known to belong to the 
PC universality class. In one spatial dimension the PC 
class is characterized by three critical exponents 
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f3 = 0.92(2), = 3.22(6), j/^ = 1.83(3) . (9) 

The two scaling exponents z^n and i>±^ are associated with 
temporal and spatial correlation lengths 



C± ^ \Pc~p[ 



(10) 



which diverge at the transition. The order parameter 
of the PC transition is the density of A-particles, corre- 
sponding to the density no of exposed sites at the bottom 
layer. Thus, in the smooth phase hq should scale as 



no {pc - p) 



(11) 



where /3o = /3 is the density exponent of the PC class. 
More generally, we expect no to obey the scaling form 



(12) 



where $o is a universal scaling function. As will be shown 
below, numerical simulations confirm the validity of this 
scaling form. This may be surprising since the ^-particles 
do not evolve independently, but are coupled to dynamic 
processes at higher levels of the interface. However, in the 
present models this feedback from higher levels is com- 
paratively weak and does not seem to affect the critical 
behavior of the ^-particles at the roughening transition. 

In order to investigate the critical behavior at higher 
levels, let us introduce the densities of sites i with hi < k 



nk 



k 



fc = l,2,. 



(13) 



Following the ideas of Ref. Q we expect ni,n2,... to 
scale in the same way as no with certain critical expo- 
nents Pk, and i^ii.fc. However, since the correlation 
lengths of neighboring layers should coincide, the expo- 
nents v±^k and J^ii.fc should not depend on k, leading to 
the scaling form 



nfe(i,L,e)^t-^^/''ii <^k(te''\ t/L') 



(14) 



Our numerical results support the validity of this scaling 
form for fc > 1 in an intermediate scaling regime. 



B. Estimation of pc and Sk 

In order to determine the critical exponents accurately, 
we perform time-dependent simulations of variant C on a 
parallel computer. Using a large lattice of 24000 sites (as- 
sociating one processor with each lattice site) we measure 
no, . . . , n2 as functions of time up to 2 • 10^ time steps av- 
eraged over 10'' samples. At criticality, the densities are 
expected to decay as 




FIG. 8. Effective exponents 5k{t) measured in model 
C with 24000 sites averaged over 10* samples. The first 
panel shows 5o{t) with p ranging from 0.3406 to 0.3409 from 
top to bottom, leading to the estimates pc = 0.3407(1) and 
So = 0.275(10). The curves for Si{t) and S2{t) at p = pc are 
shown in panel 2 and 3. The last panel shows the magnetiza- 
tion M as a function of time. 



Ukit) t 



(15) 



where Sk — I3k/v\\- In Fig. || the effective exponent 5o{t) 
(i.e., the local slope of no(i) in a log-log plot, see [^) is 
plotted against i in a logarithmic scale for various values 
of p. As can be seen, the curve that appears to become 
approximately horizontal for large t corresponds to the 
critical point pc — 0.3407(1). For the exponent i5o we 
obtain the estimate 



6o = 0.275(10) , 



(16) 



which is in agreement with the expected value of the PC 
class 5 = l3/v\\ = 0.285(5). 

The estimation of the exponents 5k at higher levels is 
more difficult. Averaging 5k{t) over two time decades we 
obtain the estimates 



^1 = 0.205(15) , (52 = 0.13(2) 



(17) 



Less accurate but compatible results for the other vari- 
ants are listed in Table |. Thus, the critical behavior at 
the first few layers is universal and does not depend on 
the RSOS constraint and the type of updates. However, 
as can be seen in Fig. |8[ the effective exponents i5i {t) and 
§2 {t) do not saturate at a constant value in the long-time 
limit. Instead they continue to increase on a logarithmic 
time scale. At present it is not clear whether the effec- 
tive exponents will eventually saturate after very long 
time. The drift of the exponents may indicate violations 
of scaling in the long-time limit which are neither related 
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to finite-size effects nor to numerical errors in the value 
of Pc- Thus, the estimates have to be taken with 
care. Similar violations of scaling have been observed in 
the monomer case and seem to be an intrinsic property 
of roughening transitions driven by absorbing-state tran- 
sitions. We will come back to this problem in Sect. [vj . 



C. Spontaneous symmetry breaking 

The growth processes considered in this paper are 
translationally invariant in height direction. This sym- 
metry is spontaneously broken in the smooth phase where 
the models select one of the heights as the bottom layer of 
the interface. In order to quantify this symmetry break- 
ing we consider an order parameter 



(18) 
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FIG. 9. Finite size results for the parallel model for 
L = 64, 128.. .4096. The density at the bottom layer no(t) 
is plotted as a function of time. The inset shows a data col- 
lapse which is explained in the text. 



where j runs over all lattice sites. This order param- 
eter can take positive or negative values and thus re- 
minds of a magnetization. In the stationary state of the 
smooth phase |M| is positive and vanishes at the rough- 
ening transition. As in the monomer case, we expect the 
magnetization to scale near the roughening transition as 



\M\ 



(19) 



where is a critical exponent. Thus, starting from a 
flat interface, \M\ should decrease at criticality as t"^/"^!! . 
Performing numerical simulations (see Fig. || and Table |) 
we find — 0.76 for both the restricted and the unre- 
stricted variants, suggesting that is a universal critical 
exponent. As will be shown in the following Section, 
this exponent can also be seen in unidirectionally cou- 
pled PC processes. Its numerical value 6 ~ 2.5 is much 
larger than in the monomer case where the value was 
found to be 0.65 Q]. At present it is not known whether 
6 is independent or related to the other bulk exponents 



D. Finite-size simulations at criticality 

The dynamic exponent z = i'wIvjl can easily be veri- 
fied b y fi nite-size simulations. According to the scaling 
form (p^, several measurements of no(i) for different sys- 
tem sizes should collapse onto a single curve if no is 
plotted versus t/L^. As shown in Fig. ^, the best data 
collapse is obtained for 5^ — 0.27 and z ~ 1.7, which is in 
agreement with the known PC exponents 5 = 0.285(5), 
z — 1.76(3), and the previous estimate of z in Sect. [V. 



E. Off-critical simulations 

In order to determine the exponent (3q directly, we per- 
form off-critical simulations in the steady state of the 
smooth phase. Here we expect the stationary density rig 
to scale as {pc—pY°, where Pq ~ 0.92(2) is the order pa- 
rameter exponent of the PC class. Surprisingly, a rough 
estimate of /3o in standard Monte-Carlo simulations of 
variants A and C yields much smaller values of about 
0.6 which deviate from the expected value by more than 
30%. Similar deviations are observed if vi^ and are 
determined in off-critical simulations. For example, the 
estimate v\\ — 2.4(5) in Sect. IV is much smaller than the 
expected PC value v\\ ~ 3.2. 

Off-critical steady state simulations are known to be 
quite inaccurate because of long transients. Moreover, 
they are extremely sensitive to errors in the estimation 
of Pc- However, in the present case the deviations have 
a different origin. This can be shown by determining 
the effective exponent /3o(e) which is defined as the local 
slope of no(e) in a log-log representation between the data 
points (i — 1, i) 



Inei - lnei_i 



where ti — Pc ~ Pi, providing an estimate 



/3o = lim/3o(e) ■ 

€^0 



(20) 



(21) 



As shown in Fig. |l^, the effective exponent increases 
slowly and tends towards the expected PC value as e ^ 0. 
A simple linear extrapolation yields the estimate 0.9 
which is in rough agreement with the expected value 
of /3o. 
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FIG. 10. Left: The effective exponent /3o(e) as a func- 
tion of e in model C with L = 24000 sites. After an equi- 
libration of 4 X 10^ time steps the bottom layer density is 
averaged over 10^ samples. A simple linear fit gives the esti- 
mate Po ~ 0.89. Right: Corresponding simulation results for 
the truncated model. Here the effective exponent is almost 
independent of e. Notice the different ranges of e. 



The drift of /3o(e) is unusual and may be related to 
the influence of dynamic processes at higher levels of the 
interface. To support this hypothesis, we analyze the 
truncated version of model C where the heights are re- 
stricted to take the values and 1 in order to eliminate 
the feedback from higher levels. Repeating the same type 
of analysis at the critical point Pc = 0.583(1) we observe 
/3o(e) to be in agreement with the expected PC value 
over a wide range of e. Thus we are led to the conclusion 
that dynamic processes at higher levels are responsible 
for corrections to scaling in the off-critical regime. 



F. Unusual scaling for random initial conditions 

So far we have considered the temporal evolution start- 
ing from an entirely flat interface. It would be natural to 
expect the asymptotic scaling behavior to be universal 
for any initial interface configuration with finite width 
and short-range correlations. Surprisingly this is not 
true. For example, starting with random initial condi- 
tions hi — 0,1 the densities rik turn out to decay much 
slower. For restricted variants we observe an algebraic 
decay of no with an exponent 5o ~ 0.13 which differs sig- 
nificantly from the value 0.275 for flat initial conditions. 
Similarly, the critical properties of the faceting transition 
ai p — 0.5 are affected by random initial conditions. 

The nonuniversal behavior for random initial condi- 
tions is related to an additional parity conservation law. 
In fact, the dynamic rules (|l])-(^ not only conserve par- 
ity of the particle number but also conserve parity of the 
droplet size. Starting with a flat interface the lateral 
size of droplets is always even, allowing them to evapo- 
rate entirely. However, for a random initial configuration, 
droplets of odd size may be formed which have to recom- 
bine in pairs before they can evaporate, slowing down the 
dynamics of the system. 

In the language of BAW2's the additional parity con- 



servation law is due to the absence of nearest-neighbor 
diffusion. Particles can only move by a combination of 
offspring production and annihilation, i.e., by steps of 
two lattice sites. Therefore, particles at even and odd 
lattice sites have to be distinguished. Only particles of 
different parity can annihilate. Starting with a fully occu- 
pied lattice all particles have alternating parity through- 
out the whole temporal evolution, leading to the usual 
critical behavior at the PC transition. For random ini- 
tial conditions, however, particles of equal parity cannot 
annihilate, slowing down the decay of the particle den- 
sity. BAW2's and other particle processes with parity- 
conserving diffusion will be analyzed in a future study. 



G. Comparison with other models 

Recently, Park and Khang proposed a growth model 
which involves two symmetric particle species p6| . Par- 
ticles are adsorbed (desorbed) with probabihty q (1 — <;) 
except at those sites where the nearest neighbors are oc- 
cupied by particles of the other (same) species. Thus, the 
dynamic rules mimic a PC transition with two symmet- 
ric absorbing states [0. However, the critical exponents 
at the roughening transition were found to differ signif- 
icantly from the expected PC values. In fact, once a 
new layer is completed, the configuration of the previous 
layer becomes frozen. Thus the kinks between different 
domains do not act as pinning centers which destroys 
the PC transition. Consequently, the interface starts to 
propagate for q > Qc- 

Another parity-conserving model introduced by Noh 
et. al. |27]] was inspired by interacting monomer-dimer 
models [|16| exhibiting a PC transition. The dynamic 
rules involve diffusion, order-preserving branching, and 
order-breaking branching. In absence of order-breaking 
branching, the interface evolves within a monolayer and 
undergoes a pre-roughening transition belonging to the 
PC class. In presence of order-breaking processes several 
layers may be formed. However, in this case the interface 
is always rough and the transition is lost. 

Even more recently, the same authors investigated a 
class of parity-conserving growth models p9| ] , generaliz- 
ing the rules (|l|)-(§) by the addition of a process for the 
evaporation of dimers from the middle of plateaus. A 
similar generalization of monomer models has been stud- 
ied in the context of nonequilibrium wetting . In both 
cases the additional evaporation process destroys the sta- 
bility of the smooth phase. In the dimer models, however, 
the roughening transition of the dimer models is replaced 
by a faceting transition between a rough and an (oppo- 



sitely) faceted phase. Thus, the models of Ref. |19| link 
anomalous roughening and faceting transitions. 
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FIG. 11. Extended particle interpretation. Dimers are 
adsorbed {2A — > 0) and desorbed {A — > 3A) at the bottom 
layer. Similar processes take place at higher levels. 



VI. UNIDIRECTIONALLY COUPLED 
BRANCHING-ANNIHILATING RANDOM 
WALKS WITH EVEN NUMBER OF OFFSPRING 




FIG. 12. Unidirectionally coupled PC processes. Left: 
Densities nA,nB,nc as functions of time. Right: Magnetiza- 
tion M{t) as a function of time averaged over 1000 indepen- 
dent reahzations. 



In order to understand the universal properties of 
parity-conserving roughening transitions let us extend 
the particle interpretation. Assuming that h = is 
the bottom layer of the interface, we associate parti- 
cles A, B,C, . . . with sites at height hi < 0, 1, 2, . . . (see 
Fig. ^l]). As explained before, the dynamic processes at 
the bottom layer may be interpreted as an effective re- 
action A 3^,2^4 0. Similarly, the B-particles on 
top of the first layer react by B ^ 3B, 2B 0. Clearly, 
the temporal evolution of the i?-particles strongly de- 
pends on the actual configuration of the A-particles. On 
the one hand, an A-particle implies the presence of a 
B-particle at the same site, giving rise to an effective re- 
action A ^ A + B. On the other hand, both the RSOS 
constraint and the restriction that dimers react at ad- 
jacent sites of equal height may inhibit the above reac- 
tions, introducing an effective feedback from higher levels 
to lower ones. However, as suggested by our numerical 
results, this type of inhibiting feedback does not affect 
the critical behavior at the roughening transition. 

Similarly, the C-particles perform an effective BAW2 
on top of the second layer. Therefore, the critical be- 
havior of the first few layers may be described in terms 
of a simplified particle model where several BAW2 pro- 
cesses are unidirectionally coupled according to the reac- 
tion scheme 

A^3A B ^3B C 3C 

2A^0 25 ^ 2C ^ (22) 

A-^A + B B-^B + C C ->C + D,... 

generalizing the concept of 'coupled DP' We propose 
this reaction scheme to characterize the universal behav- 
ior of the dimer models at the roughening transition. 

To support this hypothesis, we study the reactions ( |2^ ) 
by Monte-Carlo simulations. To this end we use three 
copies of the BAW2 introduced in Ref. Q coupled by the 
rule that a selected ^(i?)-particle instantaneously creates 
a B(C)-particle at the same position, provided that the 
target site is empty. Using a lattice size of 2500 sites we 
measure the particle densities nA , ns , nc as functions of 
time at criticality. As shown in Fig. |l^, they display 



essentially the same behavior as the densities no,ni,n2 
in the dimer model. Averaging over one decade in time 
we obtain the exponents 

Sa ^ 0.280(5) , 6b = 0.190(7) , 6c = 0.120(10) , (23) 

which are in fair agreement with the corresponding expo- 
nents of the dimer model. We also measured the analog 
of the magnetization M = ua — tib ~'r uq — njy + . . . 
in 20 coupled PC processes, obtaining an algebraic de- 
cay - t^^/^ii with the exponent e/iy^i = 0.74(3). 
This confirms the hypothesis that the roughening transi- 
tion of the dimer models belongs to the universality class 
of unidirectionally coupled PC processes. Notice that 
the coupling A ^ A + B violates parity conservation at 
the i?-level. However, using a parity-conserving coupling 
such as A ^ A + 2B we obtain similar results at the 
transition. 

Before turning to the field-theoretic formulation of cou- 
pled PC processes, let us consider the mean field approx- 
imation of the reaction scheme ( |2^ ) 

dtUA = cm A ~ Xua , 

dtUB = cms - Xn% + M'T'A , (24) 
dtnc = cmc — Xuq + iiub , ■ ■ ■ 

where nA,nB,nc correspond to the densities no,ni,n2 
in the growth models, a and A are the rates for off- 
spring production (dimer evaporation) and pair annihi- 
lation (dimer deposition), respectively. The coefRcient fi 
is an effective coupling constant between different parti- 
cle species. Since these equations are coupled in only one 
direction, they can be solved by iteration. Obviously, the 
mean- field critical point is Cc = 0. For small values of 
a the stationary particle densities in the active state are 
given by 

/ \ 1/2 / \ 1/4 

corresponding to the mean field critical exponents 



9 



/3. 



MF 
A 



1/2, 



oMF 

Pc 



1/4,. 



(26) 



These exponents should be vahd for d > dc = 2 (see be- 
low). Solving the asymptotic temporal behavior we find 
I'll = 1, implying that S^^^ — 2~^. Although this simple 
mean-field calculation does not include parity conserva- 
tion, it explains the reduced values of the exponents at 
higher levels. The different numerical values in one di- 
mension are due to fluctuation corrections which may be 
computed within a field-theoretic renormalization group 
approach. A field theory for a single BAW2 was intro- 
duced and studied in detail by Cardy and Tauber [[lit . 
Following their notation, the effective (unshifted) action 
of unidirectionally coupled BAW2's should be given by 

S[lpO,1pl,1p2, ■ ■ . ,■00,-01,-02, . . .] = 

OO 

d^xdt Y,{Mdt ~ ^V2)^fc - A(l - i>l)^l + 

T.2 



+a{l - iil)%pki'k + m(1 - '(/'fe)i/'fe-i-0fc-i| 



(27) 



n.(t) 0.1 r 




FIG. 13. Parity-conserving PNG model. The densities 
no ■ ■ .ns are shown as functions of time. The exponents 5k 
are estimated in the interval 10"^ < t < 10^, suggesting that 
the model belongs to the universality class of unidirectionally 
coupled PC processes. The inset shows the magnetization 
M{t). The dashed line indicates the slope —0.77. 



where = = 0. Here the fields tpk and ipk rep- 
resent the configurations of the system at level k. How- 
ever, even for a single copy in 1+1 dimensions the field- 
theoretic treatment poses considerable difficulties. They 
stem from the presence of two critical dimensions: d^ — 2, 
above which mean-field theory applies, and ~ 4/3, 
where for d > d'^ {d < d'^) the branching process is rele- 
vant (irrelevant) at the annihilation fixed point. There- 
fore the physically interesting spatial dimension d = 1 
cannot be accessed by a controlled e-expansion down 
from upper critical dimension dc = 2. A one-loop calcu- 
lation for the action ( p7| ) would be an even more difficult 
task. 

A fundamental problem of the field theory for coupled 
DP describing the monomer case is the relevant coupling 
/i which grows under renormalization group transforma- 
tions Q. Moreover, even in a one- loop calculation cer- 
tain infrared divergent diagrams are encountered which 
are proportional to /i. Similar difficulties are expected 
for unidirectionally coupled PC processes. They might 
be responsible for violations of scaling in the long time 
limit appearing as a curvature of and nc in Fig. |lj. 



VII. PARITY-CONSERVING POLYNUCLEAR 
GROWTH MODELS 

Roughening transitions related to coupled DP were 
first observed in so-called polynuclear growth (PNG) pro- 
cesses [^|J2^-Q. In these models the interface grows by 
nucleation processes and deterministic growth of terraces. 
The use of parallel updates ensures that the maximal 
propagation velocity is 1 . Depending on the rates for nu- 
cleation and terrace growth, the models exhibit a rough- 
ening transition from a moving rough phase to a smooth 



phase propagating at maximal velocity. As shown in 
Ref. 1^ , PNG models may be viewed as upside-down ver- 
sions of monomer models in a co-moving frame. In some 
cases it is even possible to relate PNG and monomer mod- 
els exactly to each other, resolving the apparent paradox 
that the transition in PNG models requires parallel up- 
dates whereas for monomer models the type of updates 
does not play a role. 

In order to understand the DP mechanism, let us con- 
sider the PNG model introduced by Kertesz and Wolf 0| 
which is updated synchronously in two sub-steps. At 
first all up (down) steps of the interface move determin- 
istically to the left (right) over a distance of u lattice 
spacings. Then all heights are increased with probabil- 
ity p. Thus the model is unrestricted. For u — 1 this 
dynamic rule can be expressed as a single parallel up- 
date 

where m(i) = max[/ii_i(i), /ii(t), /ii_|_i(i)]. Starting from 
a flat interface hi{0) = 0, the sites at maximal height 
hi{t) = t may be considered as the active sites of a DP 
process. The nucleation process turns active into inactive 
sites with probability 1 —p while the deterministic growth 
of terraces resembles offspring production. Therefore, if 
p is large enough, the interface is smooth and propagates 
with velocity 1. Below a critical threshold pc — 0.539(1), 
however, the growth velocity is smaller than 1 and the 
interface evolves into a rough state. As shown in Ref. [||, 
this model is a realization of unidirectionally coupled DP 
processes in a co-moving frame, where the order param- 
eters Uk are given by 
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"feW = ■ (29) 

i h=0 

In the following we introduce a parity-conserving PNG 
model which belongs to the universality class of unidircc- 
tionally coupled PC processes. The model is defined on 
a one-dimensional lattice with periodic boundary condi- 
tions and evolves by sublattice-parallel updates. In the 
first half time step pairs of sites + I) with even i are 
updated. If hi{t) ^ hi+i{t), the heights are incremented 
by one step 

h,{t + 1/2) ^ h,{t) + 1 , (30) 
h,+i{t + 1/2) ^h^+iit) + l . 

If, however, the two heights are equal, they are updated 
by the probabilistic rule 

h,{t + 1/2) = h,+i(t + 1/2)^ 

( ni{t) + 1 with prob. p , , 

1 ni{t) with prob. 1 — p ' 

where m{t) = max.[hi-i{t) , hi{t) , hi+i{t) , hi+2it)]- In the 
second half time step the same update rule is applied 
to odd pairs of sites. Clearly, this model generalizes the 
PNG model of Ref. Q and conserves parity at each height 
level in a co-moving frame. The conservation law leads 
to the formation of pinning centers moving at maximal 
velocity. 

Performing Monte-Carlo simulations we observe a 
roughening transition at the critical threshold Pc = 
0.5697(3). Starting from a flat interface, we measure the 
densities nk{t) defined in Eq. ( p9| ) as functions of time. 
As shown in Fig. the temporal decay of Uk is similar 
to the one observed in the dimer model. Averaged over 
two decades in time we obtain the estimates 

<5o = 0.28(1), (5i= 0.21(2), ^2 - 0.14(2), (32) 

which are compatible with the values listed in Table |. 
The magnetization M{t) defined in Eq. ( p^ ) does not 
show a clean power law behavior but seems to approach 
an asymptotic decay i"" *"^ in agreement with previous 
findings (see inset of Fig. |l^). Therefore, we conclude 
that the roughening transition of the parity-conserving 
PNG model belongs to the universality class of unidirec- 
tionally coupled PC processes. 



VIII. CONCLUSIONS AND OUTLOOK 

In this paper we have investigated a class of parity- 
conserving growth processes in which dimers adsorb at 
sites of equal height and desorb at the edges of ter- 
races. At a critical growth rate p = Pc the mod- 
els display a roughening transition from a smooth to a 



rough phase. In order to demonstrate the robustness of 
this transition, we have studied four variants of a 1-1-1- 
dimensional parity-conserving growth process with and 
without RSOS constraint, using either random-sequential 
or sublattice-parallel dynamics. In addition, we have in- 
troduced a parity-conserving polynuclear growth process 
where a similar transition takes place in a co-moving 
frame. In all cases the roughening transition is char- 
acterized by the same type of critical behavior. 

The investigated dimer models generalize previously 
studied monomer models. Their essential feature is a 
parity conservation law at each height level, changing the 
universal properties of the roughening transition. The 
conservation law leads to the formation of pinning cen- 
ters separating regions of even and odd parity. Thus, 
in contrast to the monomer case, the interface remains 
pinned to the initial height level. At the transition the 
width is found to increase logarithmically with time due 
to a slow diffusion and annihilation of the pinning cen- 
ters. In Eq. (0) we have proposed a finite-size scaling 
form for logarithmic roughening which is confirmed by 
high precision simulations. Moreover, we have shown 
that the universal critical behavior at the first few lay- 
ers may be described in terms of unidirectionally coupled 
branching-annihilating random walks with two offspring. 
Thus we suggest that the concept of undirectionally cou- 
pled PC processes defines a whole universality class of 
parity-conserving roughening transitions. 

In case of equal rates for dimer adsorption and desorp- 
tion all variants undergo a second phase transition where 
the width increases algebraically. In the restricted case 
this transition has been identified as a faceting transition 
from a rough to a faceted phase [|l^ . In the unrestricted 
variants, however, we observe a transition from a rough 
to a freely growing phase characterized by spiky interface 
configurations. 

There are various possible extensions and generaliza- 
tions of the models studied in this paper. Very recently, 
Noh et al. |l9| investigated a generalization of variant 
A where dimers may also evaporate from the middle of 
plateaus. Even at a very small rate this additional pro- 
cess destroys the stability of the smooth phase, turning 
it into a faceted phase. Remarkably, a sharp transition 
between the faceted and the rough phase still remains, 
leading to interesting crossover phenomena between dif- 
ferent universality classes which have not been studied 
so far. It would also be interesting to investigate parity- 
conserving growth processes in higher dimensions. Since 
the upper critical dimension d'^ is less than 2, we expect 
the roughening transition - if still existing - to be de- 
scribed by mean-field exponents. One may also consider 
growth processes of n-mers where the number of particles 
at each height level is preserved modulo n. Especially in 
higher dimensions, these n-mers might appear in differ- 
ent shapes and orientations. After all it would also be 
interesting to find experimental realizations for deposi- 
tion and evaporation of composite particles. 
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